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We study the dynamic behaviour of a quantum two-level system with periodically varying pa- 
rameters by solving the master equation for the density matrix. Two limiting cases are considered: 
multiphoton Rabi oscillations and Landau-Zener (LZ) transitions. The approach is applied to the 
description of the dynamics of superconducting qubits. In particular, the case of the interferometer- 
type charge qubit with periodically varying parameters (gate voltage or magnetic flux) is investi- 
gated. The time-averaged energy level populations are calculated as functions of the qubit 's control 
parameters. 
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I. INTRODUCTION 

The quantum two-level system is a model to describe a 
number of physical objects such as atoms, quantum dots, 
molecular magnets, etc. The system can be excited from 
the ground state to the upper state by changing in time 
its parameters (due to external fields). If the parame- 
ters vary adiabatically slow, the excitation mechanism 
is called Landau-Zener (LZ) transitionii^; if the ampli- 
tude of the field is small and its frequency is comparable 
to the level distance, Rabi oscillations occur— i^. Thus 
the occupation of the levels in a two-state system can be 
controlled by several parameters, e.g., by amplitude and 
frequency of an external field^. 

Recently, the two-level model was used to describe 
Josephson-junction systems^, i.e., both charge qubits^ 
and flux qubita^. The subject of the present work is the 
investigation of the dynamic behaviour of superconduct- 
ing qubits with periodically swept parameters, which is 
important from the point of view of state control and 
readout. We will relate our results to some other ar- 
ticles concerning the resonant excitation of a two-level 
system^" ^■^ which was shown to be relevant for supercon- 
ducting qubits tooi^~— . Particularly, we will describe 
the dynamic behaviour of the interferometer-type charge 
qubiti^-i^. 

To study theoretically the dynamic behaviour of a 
Josephson qubit, we make use of the master equation for 
the density matrix rather than the Schrodinger equation 
because this allows to take into account both relaxation 
and non-zero temperature effects. (However, in this pa- 
per we will assume the zero-temperature limit, because 
we are interested in the influence of relaxation processes 
only on the qubit dynamics.) 

The diagonal components of the density matrix define 
the probabilities of finding the system in the respective 
states of the basis in which the density matrix is pre- 
sented. Thus, executing the calculations we ought to deal 



with a particular basis. Our calculations are mostly car- 
ried out in the stationary basis {|— ) , |+)} of the eigen- 
states of the Hamiltonian in the absence of time- 
dependent terms. We do that for two reasons. First, 
it is convenient to describe Rabi oscillations and mul- 
tiphoton transitions. Second, in the case of the charge 
qubit these states {|— ) , |+)} are eigenstates of the cur- 
rent operator, which is related to the experimentally mea- 
surable valuesSi. The important point is that we can 
get the occupation probability of any state provided we 
know the probabilities for the states in a particular basis. 
To demonstrate this, we will change over from the sta- 
tionary basis to the so-called adiabatic basis (consisting 
of the instantaneous eigenstates of the time-dependent 
Hamiltonian H) to describe the LZ effect. We emphasize 
that the presented results are valid for the description of 
any two-level system with periodically swept parameters, 
particularly, of a superconducting qubit. In a Josephson- 
junction qubit the gate voltage or the magnetic flux can 
be modulated periodically. 

The paper is organized as following. In Sec. II the ba- 
sic equations are presented. In Sec. HI we study multi- 
photon processes and LZ transitions in a two-level system 
with time-dependent parameters. We apply the respec- 
tive results to the phase-biased charge qubit in Sec. IV. 



II. THE BASIC EQUATIONS 

We start from the Hamiltonian of a two-level system 
(see e.g. Ref. 6) 
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in the basis of "physical" states {|0) , |1)}, where az |0) = 
|0), dz |1) — — |1), ^x,y,z are the Pauh matrices. The 
" physical" states are the eigenstates of the Hamiltonian 



2 



for b'^x^/B^z^ 0. In the case of a charge qubit, 
these states correspond to a definite number of Cooper 
pairs on the island. For a flux qubit, they correspond 
to a definite direction of the current circulating in the 
ring. The Hamiltonian (I) is diagonalized by means of 
the matrix 



5 = cxp ( i-dy ] = 
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The eigenstates |— ) and |+) of the time-independent 
Hamiltonian H'^^'^ are connected with the initial basis: 
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Next, we introduce the time-dependent terms into the 
Hamiltonian (1), 

We consider two situations, 

(a): = Bz=BM = Bf^+B^^\t), (2) 

(6): (3) 

where the time-independent/dependent terms are 
marked with (0)/(l) indices. Making use of the trans- 
formation = S^^HS, we get the Hamiltonian in 
the energy representation {|— ) , |+)} corresponding to 
these cases: 
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AE^ BP(t) , 

- — — (sm Xtx + cos Xtz ) , (4) 



——Tz — (cos Xtx - sm Xtz) . (5) 



For convenience, we use different notations for the Pauli 
matrices, Gi and ri, which operate in the bases {|0) , |1)} 
and {|— ) , respectively. 

We emphasize that after the substitution 



B^^\t)^B(^\t), 

^(0) 5(0) 
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the problem (b) coincides with problem (a). This tran- 
sition from Eq. Q to Eq. Q corresponds to a 7r/2 
rotation about the y-axis. 



To unify expressions l@J and (O to get the equations for 
numerical calculations, we write down the Hamiltonian 
H/ as following: 



H' = —Tx H T, 

2 2 
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The quantum dynamics of our two-level system can 
be characterized within the standard density-matrix 
approach™. The time-evolution of the total system com- 
posed of the two-level system and the reservoir is de- 
scribed by the Liouville equation. After tracing over the 
reservoir variables, the Liouville equation can be sim- 
plified to the so-called master equation for the reduced 
density matrix p. It can be written in the form 
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l + Z X 
X + iY 1 
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which ensures the condition Trp = 1. The effect of re- 
laxation processes in the system due to the weak cou- 
pling to the reservoir can be phenomenologically de- 
scribed with the dephasing rate F^ and the relaxation 
rate Treiax(~)- Then the master equation takes the form 
of Bloch equationsS^: 



= -CY - T^X, 



dX 

~dt 
dY 

-J- = -AZ + CX- F^r, 
dt 

— =AY- Trela. {Z - Z{Q)) . 

dt 



(8) 
(9) 
(10) 



(Throughout the paper we admit h — 1.) From these 
equations we get Z{t) which defines the occupation prob- 
abihty of the upper level |-|-), P+{t) — P22{t) = 5(1 ^ 
Z(t)). We choose the initial condition to be X(Q) = 
Y{0) = 0, Z{Q) = 1, that corresponds to the system in 
the ground state |— ). We have calculated the time evo- 
lution of P+{t) (which is essential e.g. for the snap-shot 
measurements^^) as well as the time-averaged probabil- 
ity P+ (which is essential e.g. for the impedance mea- 
surement techniques^). We note that the asymptotic ex- 
pression, -P+(i)|(^QQ, is periodic in time with the period 
T = 27r/cj of the time-dependent term of the Hamiltonian 
H^^\t) (see Sec. 12 of Ref. 5). 



III. NON-STATIONARY EFFECTS IN A 
TWO-LEVEL SYSTEM 

A. Rabi oscillations 

Hereafter we will treat the problem (a) (see Eq. Q), 
making use of the following notations: 

H = Adx + x{t)az, x{t) — Xof f + XQshiijjt (11) 
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This reformulation is convenient for comparing our re- 
sults with the results of other papers^"— . Then the 
Hamiltonian can be rewritten: 



H' = 2~^^ ^ sinwt • V, 

V - _ ^^°// - 

^ ~ AS^" AE 



(13) 



x„M=10, ci)/A=0.1 



(b) 
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First consider the situation ^o// = 0, then the differ- 
ence between the stationary energy levels is /S.E ~ 2A. 
In the case 



Ac 



AE -^AE and xq < AE, (14) 



one can use the so-called rotating wave approximation, 
and the result is'^^ 

P.it) - ^^^^ (1 - cos ^^IHA.)h) . (15) 

For the average probability, it follows 



P+ = 



1 xl 



2x1 + i^^f 



(16) 



This means that at a; = A_E there is resonance, Pj. 



+ - 5' 

and P+{t) is an oscillating function with the frequency 
Xq. This is illustrated in Fig. The width of the peak 
at cij = AE of the P^-lo curve at the half-maximum (i.e., 
at P+ — 1 /4) is approximately 2xq (see the upper panels 
of Fig. El. 

Resonant excitations of a two-level system for 
xq/AE <C 1 may occur not only at w = AE. In the 
K-th order approximation, {xq/ AE)^ , resonances exist 
at w ~ AE / K'^^'^T-. For Xoff = 0, the resonances occur at 
odd K only. The dependence of P+ on time in the vicin- 
ity of K-th resonance is described by a relation similar to 
Eq. (Uni) (with substitutions Auj Auj^'^'^ ~ Kuj - AE 
and Xq — > Xq). Hence, the width of the K-th resonance 
is of the order of x^ . In the resonance, Aw^^' = 0, the 
energy levels |— ) and |+) are equally populated so that 



P+iu 



OJ 



AE/K) = 1/2, which is the maximum 



value of P+ , since the population inversion in a two-level 
system is not possible. 

When Xoff 9^ 0, the matrix V has non-zero diagonal 
elements (see Eq. ^^). The appearance of the non-zero 
diagonal elements does not influence Rabi resonances at 
odd K , but rather results in the generation of the reso- 
nances at even K. This is demonstrated in Fig. |21 and 
discussed in Sec. III.C. 



FIG. 1: Time dependence of upper-level occupation probabil- 
ities, (a) Rabi oscillations in P+ with the period Tr = 2-k/xo, 
(b) LZ transition in Pf (see Sec. III.B), (c) and (d) P+ prob- 
ability evolution in the case of periodically swept parameters 
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FIG. 2: Dependence of the probability P+ on the frequency 
u) for different xq at r0 = T relax = and at Xoff = (solid 
line) and Xoff = 0.2 A (dashed line). (Only the first few 
resonant peaks are plotted; the others, which are very narrow, 
are not shown at the graphs.) Inset: enlargement of the low- 
frequency region. 



B. LZ effect 

The LZ effect manifests in the non-adiabatic transition 
with the probability 



Plz = cxp - 



7rA2 

LUXq 



(17) 



between two adiabatic energy levels during a single-sweep 
evenlji. 

The time-dependent Hamiltonian H is diagonalized in 
the adiabatic basis, denoted as {||) , |t)}, by the 5*1 ma- 
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trix, 
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— sin cos ^ 



smr; 
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2x(t) 
' AEi '■ 



AEi{t) = 2y^A^ + xity 



The instantaneous eigenvalues of H are — 
±AEi{t)/2. 

We now can change over from the stationary basis to 
the adiabatic one, 
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Assuming Xoff — 0, we obtain 
1 
71 



S2 



cos ^ + sm ^ 

nr\Q IL — oin H 



■ cos ■ 



sinf 



Thus, provided we calculate the density matrix in the 
stationary basis p, we find it in the adiabatic basis Padiab, 



Padiab 



S2 ^pS2. 



(18) 



The initial condition for p can be obtained from the initial 
condition for Padiab by inverting the relation 1)18(1 . 

Let us now consider as an illustrative example the one- 
sweep process, 



t e 



4' 4 j ' ~' 



x{t) — xq sin ujt, x{t) e {—Xq,Xq)^ 
xo > A > tj, 

which corresponds to the LZ model—. We choose the 
initial condition to be 
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which means that the system at i = — T/4 is in the lower 
adiabatic state ||). We look for the occupation probabil- 
ity P| — Pa^iabi^) of upper adiabatic level |t), which 
equals to the LZ probability at the end of the sweep^^, 
P|(T/4) = Plz- Thus, for the functions X,Y,Z, intro- 
duced in Sec. II, we get the initial condition 



Xo 
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and find 
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This is illustrated in Fig. QJd which is equivalent to Fig. 
3d in Ref. 10. 

In the general case, Xo// ^ 0, we calculate S2 and then 
padiab according to Eq. (|18|l . Then the probability is 
given by 

P,it) . i {^'^< t)x^ 

2 2v/A2 + 2;(t)2 /a2 ' -2 



"off 



{Ax{t) - Axo 
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2y/A^+x(t)\ A 



i2 I ^2 



-.X{t). 
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This should be supplemented with the respective initial 
condition. 



C. Crossover from multiphoton Rabi resonances to 
LZ-interferometry: numerical results 

Now making use of the numerical solution of Eqs. (jSj- 
llOfl for the Hamiltonian l|13(l . we study the dependence 
of the time-averaged probability P_(. on frequency oj and 
amplitude Xq. For small amplitudes, xq ^ AE, there are 
resonant peaks in the P^-lo dependence at a; ~ AE/K, 
as described in Sec. III. A and illustrated in Fig. |21 
With increasing amplitude xq, the resonances shift to 
higher frequencies. For Xoff = 0, the resonances ap- 
pear at "odd" frequencies {K ~ 1,3,5,...) only, as it 
was studied in Ref. 9. For Xoff ^ there are also res- 
onances at "even" frequencies {K = 2,4,...), which is 
demonstrated in Fig. [3 We note that Fig. |21 is plot- 
ted for the ideal case of the absence of decoherence and 
relaxation, = Treiax = 0, when the resonant value 
is P+ — 1/2. The effect of finite dephasing, F^ 7^ 0, 
and relaxation, Treiax 7^ 0, is to decrease the resonant 
values of P+ and to widen the peaks for F^ > Treiax- 
Thus from the comparison of the theoretically calculated 
resonant peaks with the experimentally observed ones, 
the dephasing F^ and the relaxation rates Treiax can be 
obtainedi^. 

When the system is driven slowly, lo ^ A, and with 
large amplitude, xq S> A, the LZ excitation mecha- 
nism is relevant for the description of the system dy- 
namics. In the previous subsection we considered the 
LZ transition for a single-sweep event. Now we study 
the periodical driving of the system: Interferences be- 
tween multiple LZ transitions happen, which leads to 
resonant excitationsHii^iS. We will compare these res- 
onances with the multiphoton Rabi ones. 

First, we note that the resonance positions depend on 
the amplitude; we demonstrate that for lo ^ A < xq 
in the inset of Fig. |21 But it is more illustrative to 
study the resonance properties via the dependence of P4. 
on the ratio xq/lu (or, more precisely, on the phase the 
state vector picks up per half-period^2 — 4^;^ /i^) ^ jt jg 
reasonable to carry out the calculation for a fixed value 
of the product xqlo^ which in its turn defines the LZ- 
probability (see Eq. H17|l 'l. This allows us not only to 
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FIG. 3: Dependence of the probability P+ on (j) = Axq/uj at 
xqw/A^ = 0.45 (which corresponds to Plz = 10~^); = 
relax = and Xoff =0 for graphs (a) and (b); F^/A = 
Freiax/A — and Xoff//^ = 0.05 for graphs (c) and (d). 



FIG. 4: The dependence of the probabihty P| on the phase 



IV. MULTIPHOTON EXCITATIONS IN THE 
INTERFEROMETER- TYPE CHARGE QUBIT 



compare our results with the results of Ref. 12, where the 
periodicity of the resonances in (j) was predicted, but also 
to demonstrate the transition from the multiphoton Rabi 
resonances to the ones induced due to the interference of 
LZ excitations by means of the P+-4> dependence. 

In Fig. inia,b) we plotted the dependence of P+ on 
<f> for Xoff = and — T relax = 0. The first few 
peaks, in the region xq < A, corresponding to the mul- 
tiphoton Rabi resonances are situated at (/) ~ K'^ (that 
follows from the resonance relation lo — AE/K). With 
increasing (j) (which is proportional to xq) we observe the 
overall rising of the curve in accordance with the conclu- 
sions of Ref. 9. At (j) ^ 1, i.e., at iv <^ A <^ xq, the 
resonance position is 27r-periodic in agreement with Ref. 
12. The non-zero offset Xoff 7^ results in the appear- 
ance of additional resonances between the basic ones^^, 
which is demonstrated in Fig. |SIc,d). Such a feature is 
similar to the multiphoton Rabi resonances. Nonzero de- 
coherence decreases and widens the peaks we have also 
demonstrated in Fig. |3Ic,d). 

In Fig. ^we present the dependence of P| on (f> for the 
calculations carried out in the adiabatic basis (cf. Sec. 
III.B). In Fig. 01a,b) we plotted such a dependence for 
different values of at Xoff = and F^ = Treiax — 0. For 
(j)^ 1 (see Fig. 0J)) the resonances are 27r-periodic, and 
the peaks are situated at the integer values of (/>/27r, as 
predicted in Ref. 12. We also plotted the dependence of 
P| on (j) for the non-zero offset Xoff ^ (see Fig. Hc,d)), 
namely for 0o// = 7r/2. For this value the additional 
peaks become as high as the basic ones. 



A. Interferometer-type charge qubit 

In this section, we consider the quantum dynamics of 
the interferometer-type charge qubit with periodically 
varying control parametersiS"— . This qubit consists 
of two Josephson junctions closed by a superconduct- 
ing ring. The charge en of the island between the junc- 
tions is controlled by the gate voltage Vg via the capaci- 
tance Cg. The junctions are characterized by Josephson 
energies Eji, Ej2 and phase differences ipi, ip2- The 
relevant energy values are the island's Coulomb energy, 
Ec — e^/2Ctot, where Ctot is the total capacitance of the 
island, and the effective Josephson energy 

ej = {Eji + E% + 2E.nEj2 cos S) . 

An important feature of the qubit is that its Josephson 
energy is controlled by the external magnetic flux $e 
piercing the ring. In this paper, the ring inductance L 
is assumed to be small. Then the total phase difference, 
S — ipi + ip2, approximately equals to 6e ~ 27r$e/*I'oj ^^nd 
thus ej = ej{6) ~ ej{Se)- 

Within the two-level model with the basis of "charge 
states" |0) and |1) corresponding to the number of 
Cooper pairs on the island, the Hamiltonian of the 
interferometer-type charge qubit can be written aam 

H = -]^e.jd^-]^EchSz. (19) 

where Ech — ^Ec{l- — Ug) and Ug — CgVg/e. Here the 
domination of the Coulomb energy of a Cooper pair AEc 
over the coupHng energy ej is assumed, AEc/sj > 1. 
The eigenstates, {|— ) , of this Hamiltonian are dis- 
criminated by the direction of the supercurrent in the 
21 

rmgai. 



6 



The qubit is considered to be coupled not only to the 
gate but also to the tank circuit that enables both phase 
control and readout. Thus there are two possibilities to 
make the Hamiltonian of the two-level system (|19|) time- 
dependent. First, the gate voltage Vg can be driven, 



"■a 



n)p sin wi, 



(20) 



and, second, the DC and AC components of the current 
in the tank circuit can induce a periodically varying mag- 
netic flux. 



5 — 5dc + 5 AC sa\ujt. 



(21) 



Here we restrict our consideration to the case of the si- 
nusoidal time dependence of the parameters. 

Further we study the time-averaged occupation proba- 
bility of the excited state P+ . We note that because the 
two states |— ) and |-|-) belong to oppositely circulating 
currents in the ring, they correspond to different signs 
of the qubit's Josephson inductance. The latter can be 
probed by the impedance measuring techniqueS^, which 
makes it possible to observe the resonant behaviour of 
P-(- studied in the following subsection. 



B. Resonant excitation of the interferometer- type 
charge qubit 

The case of the excitation of the interferometer-type 
charge qubit via the gate (Eq. H2UI) ') can easily be related 
to the one considered in Sec. Ill; see the Hamiltonians 
and (|19|l with the diagonal time-dependent param- 
eters defined by Eq. H2U|I . Then the ratio xq/A is given 

by 



xo 
A 



rig . 



Thus, both mechanisms considered in Sec. Ill can be 
realized: multiphoton excitations {xq/A <C 1) and LZ 
interferometry (xq/A ^ 1). 

Let us now consider the second possibility of the exci- 
tation of the qubit by the time-dependent magnetic flux. 
The Hamiltonian of the interferometer-type charge qubit 
((T^ with the periodically varying phase 6 (Eq. (PT|l ) is 
related to the Hamiltonian of a two-level system consid- 
ered in Sec. II (Eq. ||SJ)) by the following relations: 



Namely, for 
S — TT :$> - — ^^==^, Sac < Sdc — tt <c tt 



V EjiEj2 



we have 



Bi^^t) ~ ^Ej^Ej25ac since;*. 



and for 



(5 - vr < ' ,„ „ -, Sac < Sdc - tt < tt 



it follows 



VEjiEj2 



\Eji — Ej2\ 

EjiEj2 
\Eji-Ej2\ 



{Sdc - tt)Sac sirnvt. 



(23) 

(24) 
(25) 

(26) 

(27) 
(28) 



When the relation (|22|l takes place, the present prob- 
lem (of the charge qubit with time-dependent magnetic 
control) can be related to the problem considered above 
in Sec. HI by the expressions ©. Then we can estimate 
the ratio xq/AE, which characterizes the mechanism of 
the excitation of the qubit. E.g., at xq/AE < 1 we ex- 
pect multiphoton resonances in the qubit's response to 
the external alternating magnetic flux. We consider this 
case below in detail. 

We note that the width of the resonant peaks is defined 
by the non-diagonal components of the Hamiltonian 
i.e. by the product xqA. Then looking at the Eqs. 125|) 
and H28|l . we conclude that the width of the resonances 
is defined by the product SAci^ " ^^s)- 

Finally, we illustrate the multiphoton resonant excita- 
tions of the interferometer- type charge qubit. Making use 
of the numerical solution of the master equation (Sec. II), 
we find the time-averaged probability P+ plotted in Fig. 
|31 The position of the multiphoton resonant peaks is de- 
fined by the relation uj — AE/K, where AE = AE{Sdc) 
is supposed to be fixed. Alternatively, when the Sdc 
component of the phase is changed and the frequency lo 
is fixed, a similar graph can be plotted with resonances 



at Sdc — S'^j^c defined by the relation AE{S)^(l) = Kuj. 



CONCLUSIONS 



Px(t) 


= £J, 






= B.M 


— Bx\sy^c=o 




- Bx{t) 


-s.(o). 


P(0) 


= Ech- 





In two limiting cases, the expressions can be simplified 
to result the time-dependent term in the form 



B'^^\t) cx smut. 



(22) 



We have studied the dynamic behaviour of a quantum 
two-level system subjected to periodical sweeping of its 
parameters. The energy levels population was calculated 
by solving the master equation for the density matrix 
with relaxation terms. Studying the population of the 
excited state in both stationary and adiabatic bases, we 
analyzed some features of the multiphoton Rabi and LZ 
effects. Particularly, we have shown certain similarities 
of the multiphoton resonances at xq ^ w with the reso- 
nances at S> A 3> w due to the interference between 
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FIG. 5: Dependence of the probability P+ on the frequency 
uj for the phase-biased charge qubit at Ug — 0.95, Eji/Ec ~ 
12.4, Ej2/Ec = 11, r^/Ec = 5 ■ 10"*, r„,,,/Bc = 10-^ 
(5ac = 0.27r, Sdc = tt + 0.27r. 



multiple LZ transitions. Based on the solution of the 
master equation for the density matrix, we described in 
detail the dynamic behaviour of the interferometer-type 
charge qubit subjected to periodically changing gate volt- 
age or magnetic flux. 
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